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The e + p charged-current deep inelastic scattering cross sections, da/dQ 2 for Q 2 be- 
tween 200 and 60000 GeV 2 , and da/dx and da/dy for Q 2 > 200 GeV 2 , have been mea- 
sured with the ZEUS detector at HERA. A data sample of 47.7 pb _1 , collected at a 
center-of-mass energy of 300 GeV, has been used. The cross section da/dQ 2 falls by a 
factor of about 50000 as Q increases from 280 to 30000 GeV 2 . The double differential 
cross section d 2 a jdxdQ 2 has also been measured. A comparison between the data and 
Standard Model (SM) predictions shows that contributions from antiquarks (u and c) 
and quarks (d and s) are both required by the data. The predictions of the SM give 
a good description of the full body of the data presented here. A comparison of the 
charged-current cross section da/dQ 2 with the recent ZEUS results for neutral-current 
scattering shows that the weak and electromagnetic forces have similar strengths for Q 2 
above M w , M§. A fit to the data for da/dQ 2 with the Fermi constant Gf and Mw as free 

parameters yields G F = (l.l71±0.034 (stat.)^;^ (syst.)+°;°}|? (PDF)) x lO" 5 GeV~ 2 and 

M w = 80.8^1 (stat.)ll;3 (syst.)+};3 (PDF) GeV. Results for M w , where the propagator 
effect alone or the SM constraint between Gf and Mw have been considered, are also 
presented. 
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1 Introduction 



Deep inelastic scattering (DIS) of leptons on nucleons is the key source of information for 
the development of our understanding of the structure of the nucleon. In the Standard Model 
(SM), charged- current (CC) DIS is mediated by the exchange of the W boson (see Fig. |l|(a)). In 
contrast to neutral-current (NC) interactions, where all quark and antiquark flavors contribute, 
only down-type quarks (antiquarks) and up-type antiquarks (quarks) participate at leading 
order in e + p (e~p) CC DIS reactions. Therefore, CC DIS provides a powerful tool for the 
flavor-specific investigation of parton momentum distributions. Even though CC events are 
kinematically less constrained than NC events due to the unobserved final-state neutrino, they 
can be identified with little background at HERA. 

First measurements of the CC DIS cross section at HERA, reported previously by the 
HI [1], |2| and ZEUS 0, f|] collaborations, extended the coverage of the kinematic range compared 
to that of the fixed-target neutrino-nucleus scattering experiments || by about two orders of 
magnitude in the four-momentum transfer squared (— Q 2 ). These analyses were based on e~p 
and e + p data samples of approximately 1 pb _1 and 3 pb^ 1 , respectively. The cross section at 
high Q 2 demonstrated, for the first time, the presence of a space-like propagator with a finite 
mass, consistent with that of the W boson. 

This paper presents results from ZEUS on the CC e + p DIS differential cross sections da/dQ 2 , 
da/dx, da/dy and d 2 a/dxdQ 2 for Q 2 > 200 GeV 2 , and comparisons to SM predictions. The 
measurements are based on 47.7 pb" 1 of data collected with the ZEUS detector from 1994- 1997 
during which HERA collided 27.5 GeV positrons with 820 GeV protons, yielding a center-of- 
mass energy ^/i = 300 GeV. The 16-fold increase in the luminosity compared to the previous 
measurements allows the double differential cross section d 2 a jdxdQ 2 to be determined in this 
high-Q 2 regime for the first time. A recent publication presented NC cross sections from 
the same data sample ||. These data, together with those presented here, permit a precise 
comparison of CC and NC cross sections up to Q 2 values of about 2-10 4 GeV 2 . 

2 Standard Model prediction 

The electroweak Born cross section for the reaction 



e^p 



v P X 



(1) 



can be written as 



dx dQ 2 



G 2 p 
Attx 



M 2 W 



M 2 W 



+ Q 2 



V + F 2 cc (x,Q 2 )-V_xF 3 cc (x,Q 2 ) 



J 2 FE c (x,Q 2 ) 



(2) 



where Gp is the Fermi constant, M\y is the mass of the W boson, x is the Bjorken scaling 
variable, y = Q 2 /xs and Y± = 1 ± (1 — y) 2 . The center-of-mass energy of the positron-proton 
collision is given by y/s = 2y / E e E p , where E e and E p are the positron and proton beam 
energies, respectively. The structure functions _F 2 and xF^ c , in leading-order (LO) QCD, 
measure sums and differences of quark and antiquark parton momentum distributions 0. For 
longitudinally unpolarized beams, 



^cc 



x[d(x, Q 2 ) + s(x, Q 2 ) + u(x, Q 2 ) + c(x, Q 2 )}, 



(3) 



1 



xif c = x[d(x, Q 2 ) + s(x, Q 2 ) - u(x, Q 2 ) - c(x, Q 2 )], (4) 

where d(x, Q 2 ) is, for example, the parton distribution function (PDF) which gives the number 
density of a down quark with momentum fraction x in the proton. Since the top quark mass 
is large and the off-diagonal elements of the CKM matrix are small, the contribution from the 
third generation quarks to the structure functions may be safely ignored M. The chirality 
of the CC interaction is reflected by the factors Y± multiplying the structure functions. The 
longitudinal structure function, F£ , is zero at leading order but is finite at next-to- leading- 
order (NLO) QCD. It gives a negligible contribution to the cross section except at y values 
close to 1, where it can be as large as 10%. 

The electroweak radiative corrections to (fj) receive contributions from initial state photon 
radiation, fermion and boson loops, and the exchange of multiple intermediate vector bosons. 
The effects of these radiative corrections are taken into account to leading order ||, so that the 
quoted cross sections in this paper are corrected to the electroweak Born level. Equation (§) is 
evaluated with G F = 1.16639 x 1(T 5 GeV~ 2 and M w = 80.41 GeV Q. The uncertainties in 
the electroweak parameters have a negligible effect both on the calculated cross sections and 
on the radiative corrections. 

Thus, the main uncertainty in the SM cross-section prediction comes from the PDF un- 
certainties, which are discussed in detail in || [ll| and are taken into account in the CC 



cross-section calculation. The resulting uncertainty in da/dQ 2 , for example, ranges from 4% 
at Q 2 = 200 GeV 2 to 10% at Q 2 = 10000 GeV 2 , and increases further at higher Q 2 . The large 
uncertainty at high Q 2 is due to the <i-quark density which is poorly constrained at high x by 
the experimental data. 



3 The ZEUS experiment 



ZEUS |L2] is a multipurpose magnetic detector designed to measure ep interactions at HERA. 
The primary components used for this analysis are the compensating uranium-scintillator 
calorimeter (CAL), the central tracking detector (CTD), and the luminosity detector. 

The ZEUS coordinate system is right-handed with the Z axis pointing in the direction of 
the proton beam (forward) and the X axis pointing horizontally toward the center of HERA. 
The polar angle 9 is zero in the Z direction. 

Tracking information is provided by the CTD |13 operating in a 1.43 T solenoidal magnetic 



field. The interaction vertex is measured with a typical resolution along (transverse to) the 
beam direction of 0.4 (0.1) cm. The CTD is used to reconstruct the momenta of tracks in the 
polar angle region 15° < 9 < 164°. The transverse momentum (p t ) resolution for full-length 
tracks can be parameterized as cr(p t )/p t = 0.0058 p t © 0.0065 © 0.0014/p 4 , with p t in GeV. 

The CAL covers 99.7% of the total solid angle. It is divided into three parts with 
a corresponding division in 9 as viewed from the nominal interaction point: forward (FCAL, 
2.6° < 9 < 36.7°), barrel (BCAL, 36.7° < 9 < 129.1°) and rear (RCAL, 129.1° < 9 < 176.2°). 
Each section is subdivided into towers which subtend solid angles between 0.006 and 0.04 
steradian. Each tower is longitudinally segmented into an electromagnetic (EMC) and one 
(RCAL) or two (FCAL, BCAL) hadronic sections (HAC). The electromagnetic section of each 
tower is further subdivided transversely into two (RCAL) or four (BCAL, FCAL) cells. Under 
test beam conditions the calorimeter resolutions were a/E = 18%/ a/ E(GeV) for electrons 
and a/E = 35% / y^W{GeV) for hadrons. The calorimeter has a time resolution of better 



2 



than 1 ns for energy deposits above 4.5 GeV. The position of the interaction vertex along the 
beam direction can also be reconstructed from the measured arrival time of energy deposits 
in FCAL 



15 



The resolution is about 9 cm for events with FCAL energy above 25 GeV and 
improves to about 7 cm for FCAL energy above 100 GeV. 

An instrumented-iron backing calorimeter [16j (BAC) measures energy leakage from the 
CAL. The muon chambers in the forward [0, barrel and rear [17] regions are used in this 
analysis to detect background events induced by cosmic-ray or beam-halo muons. 

The luminosity is measured using the Bethe-Heitler reaction ep — > ejrf [18]. The result- 
ing small angle energetic photons are measured by the luminosity monitor, a lead-scintillator 
calorimeter placed in the HERA tunnel 107 m from the interaction point in the positron beam 
direction. 



4 Monte Carlo simulation 

Monte Carlo simulations (MC) are used to determine the efficiency for selecting events, to 
determine the accuracy of kinematic reconstruction, to estimate the background rate and to 
extrapolate measured cross sections to the full kinematic phase space. A sufficient number of 
events is generated to ensure that errors from MC statistics can be neglected. The MC samples 
are normalized to the total integrated luminosity of the data. 

The ZEUS detector response is simulated with a program based on geant |19j . The gener- 
ated events are passed through the simulated detector, subjected to the same trigger require- 
ments as the data, and processed by the same reconstruction programs. 

The underlying distribution of the Z-coordinate of the event vertex is determined using a 
minimum-bias sample of low-Q 2 neutral-current DIS events as discussed in detail in ||. 

CC DIS events including radiative effects are simulated using the HERACLES 4.5.2 || pro- 
gram with the DJANGQ6 2.4 P0[ interface to the QCD programs. In HERACLES, corrections 



for initial-state radiation, vertex and propagator corrections, and two-boson exchange are in- 
cluded. The QCD cascade and the hadronic final state are simulated using the color-dipole 



model of ARIADNE 4.08 [21] and, as a systematic check, the meps model of lepto 6.5 [22 



Both programs use the Lund string model of jetset 7.4 |23| for the hadronization. A set of NC 
events generated with DJANGO is used to estimate the NC contamination in the CC sample. 
Photoproduction background is estimated using events simulated with herwig ||24|| . The 



background from W production is estimated using the EPVEC |25] generator, and the back- 
ground from Bethe-Heitler production of charged-lepton pairs is generated with the lpair [2C 
program. 



5 Reconstruction of kinematic variables 

The principal signature of CC DIS events at HERA is the presence of a large missing transverse 
momentum, Ifr. This is illustrated in Fig. |l|(b), where an event from the final CC DIS sample 
is shown. The struck quark gives rise to one or more jets of hadrons. The energetic final-state 
neutrino escapes detection, leaving a large imbalance in the transverse momentum observed in 
the detector. Ifx is calculated as 




3 



where the sum runs over all calorimeter energy deposits Ei (uncorrected in the trigger, but 
corrected in the offline analysis as discussed below), and Q{ and fa are their polar and azimuthal 
angles as viewed from the interaction vertex. The hadronic polar angle, 7^, is defined by 

fa 2 - S 2 

where 

5 = Y,(E l - EiCosOi) = Y,(E - Pz ) t . (7) 

i i 

In the naive Quark Parton Model, 7/, gives the angle of the struck quark. Another variable 
used in the selection is the total transverse energy, Ex, given by 

E T = ^E i wy6 i . (8) 



The kinematic variables are reconstructed using the Jacquet-Blondel method |27|. The 
estimators of y, Q 2 and x are: 

y JB = 5/(2E e ); Q 2 JB = If T 2 /(l - y JB ); x JB = Q 2 B /{sy JB ). (9) 

For the offline determination of fa, 5 and Et, methods developed and tested for the NC 
cross section determination || are used. The calorimeter cells with energy deposits are grouped 
into units called clusters. For each cluster, corrections depending on the cluster energy and 
angle are made for hadronic energy loss in inactive material in front of the calorimeter. The 
correction algorithm, which is based on MC, has been verified using the highly constrained 
NC events measured in the ZEUS detector. Energetic hadron jets in the FCAL direction may 
produce particles backscattered into the BCAL or RCAL (albedo). Also, particles may be 
redirected by the material between the interaction point and the calorimeter. Such effects, 
which create biases in the measurement of 7^, are suppressed by removing low energy clusters 
at polar angles much larger than the calculated value of 7/1. 



6 Event selection 

CC DIS candidates are selected by requiring large fa and a reconstructed event vertex consis- 
tent with an ep interaction. The main sources of background affecting the CC event selection 
are processes like NC DIS and high-i?^ photoproduction, where the finite resolution or en- 
ergy escaping detection in the CAL cause fa. Events not originating from ep collisions such 
as beam-gas interactions, beam-halo muons or cosmic rays can also cause substantial appar- 
ent imbalance in the transverse momentum and constitute other sources of background. The 
selection criteria described below are imposed to separate CC events from the background. 

The events are classified first according to 70, the value of 7^ measured with respect to the 
nominal interaction point. If 70 is sufficiently large, i.e. in the central region, tracks in the 
CTD are used to reconstruct the event vertex, which strongly suppresses non-ep backgrounds. 



The selection procedure designed to select these events is described in Sect. |T^. On the other 
hand, if 70 is small, i.e. in the forward region, the hadronic final state of such CC events is 
often outside the acceptance of the CTD, and thus calorimeter timing is used for the vertex 
reconstruction. The algorithm designed specifically to select such events, which tend to be at 
high x values, is described in Sect. |6.3| . The kinematic quantities are finally recalculated using 
the Z-coordinate of the event vertex (Zvtx) determined from either CTD tracks or calorimeter 
timing, depending on 70. 
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6.1 Trigger selection 



ZEUS has a three-level trigger system ||12|| . At the first trigger level, events are selected using 
criteria based on the energy, transverse energy and missing transverse momentum determined 
by the calorimeter . Generally, events are triggered with a lower threshold of these values in 



coincidence with at least one CTD track, while a higher threshold is necessary for events with 
no CTD track. The latter events have a hadronic final state boosted forward outside the CTD 
acceptance. Typical threshold values are 5 GeV (8 GeV) in missing transverse momentum, or 
11.5 GeV (30 GeV) in transverse energy, for events with (without) CTD tracks. 

At the second level, timing information from the calorimeter is used to reject background 
events inconsistent with the bunch-crossing time. Also, the missing transverse momentum is 
available with better resolution than at the first level, so that a tighter cut of 6 GeV (9 GeV 
without CTD track) can be made. 

At the third level, track reconstruction and vertex finding are performed and are used to 
reject candidate events with a vertex that is inconsistent with an ep interaction. The thresholds 
on the trigger quantities are lower than the cut variables used in the offline analysis. 

6.2 Offline selection based on a CTD vertex 

Events with 70 > 23° are required to contain a vertex reconstructed from CTD tracks and to 
satisfy the following criteria: 

• #r > 12 GeV and If T ' > 10 GeV 

tfr' is the missing transverse momentum calculated excluding the FCAL towers closest 
to the beam hole. The I?t' cut strongly suppresses beam-gas background events while 
maintaining high efficiency for CC events. 

• |^vtx| < 50 cm 

A vertex reconstructed by the CTD is required to be within the range consistent with the 
ep interaction region. 

• Tracking requirement 

At least one track associated with the event vertex must have transverse momentum in 
excess of 0.2 GeV and a polar angle in the range 15° to 164°. 

• Rejection of photoproduction 

Photoproduction events tend to have azimuthally symmetric hadronic energy flow. At 
high Et, a relatively small imbalance due to resolution effects or escaping particles can 
lead to non-negligible I?t- These events are rejected by the following cuts: Ifx / E T > 0.4 
is required for events with 20 < rfr < 30 GeV; I?t / Et > 0.55 is required for events 
with Ifx < 20 GeV. No Ifx / Et requirement is imposed on events with Ifr > 30 GeV. 
In addition, the difference between the direction of the (P x , P y ) vector calculated using 
CTD tracks and that obtained using the calorimeter is required to be less than 1 radian 
if tf T < 20 GeV and less than 2 radians if If T > 20 GeV. 

• Rejection of NC DIS 

NC DIS events in which the positron or jet energy is poorly measured can have a large 
]?t- To identify such events, a positron-finding algorithm which selects isolated electro- 
magnetic clusters pi is used. Candidate positron clusters within the CTD acceptance 



are required to have an energy above 4 GeV and a matching track with momentum larger 
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than 25% of the cluster energy. Clusters with 9 > 164° are required to have a transverse 
momentum exceeding 2 GeV. Events with a candidate positron satisfying the above cri- 
teria and 5 > 30 GeV are rejected; for contained NC events, 5 peaks at 2E e = 55 GeV. 
This cut is applied only for events with ]fc < 30 GeV. 

• Rejection of non-ep background 

Beam-gas events typically give calorimeter arrival times which are inconsistent with the 
bunch-crossing time. Such events are rejected. A muon-fmding algorithm based on 
calorimeter energy deposits or muon-chamber signals is used to reject events produced by 
cosmic-ray or beam-halo muons. 

6.3 Offline selection without CTD vertex 

Events with 70 < 23° are not required to have an event vertex reconstructed from CTD tracks. 
They must satisfy the following criteria: 

• #r > 14 GeV and If T ' > 12 GeV 

Relaxing the requirements on tracking and the CTD vertex results in an increase of 
non-ep background. To compensate for this, the requirements on the missing transverse 
momentum are tightened. 

• |Zvtx| < 50 cm 

Zvtx is reconstructed from the measured arrival time of energy deposits in FCAL. The 
relation between the timing measurement and Zvtx was determined using a large data 
sample of NC DIS events, in which a reliable Zvtx estimate can be obtained from the 
positron track even if the hadronic system is boosted in the very forward direction. 

• Rejection of photoproduction 

TpT I Et > 0.6 is required for events with Ifr < 30 GeV. This cut also suppresses beam-gas 
interactions. 

• Rejection of non-ep background 

The same timing and muon- rejection cuts are used as described in Sect. A class of 
background events which are especially troublesome in this selection branch arises from 
beam-halo muons interacting inside the FCAL. To reduce this background, topological 
cuts on the transverse and longitudinal shower shape are imposed; these reject events 
where the energy deposits are much more strongly collimated than for typical hadronic 
jets. Another characteristic of muons traversing the detector parallel to the beam line is 
a coincidence of energy deposits in the RCAL and FCAL at similar (X, Y) positions. If 
such a coincidence is found, the event is rejected. 

NC DIS is negligible in this selection branch. 

6.4 Final event sample 

In order to restrict the sample to regions where the resolution in the kinematic variables is 
acceptable and the background is small, further requirements Q 2 JB > 200 GeV 2 and i/jb < 0.9 
are imposed. The cross sections presented below are corrected to the full y range using the SM 
^/-dependence described by (^). 
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The combined selection efficiency of the above cuts for most of the x and y region (0.1 < x, 
0.1 < y < 0.9) is typically 90%. At low x or high y, the efficiency decreases due to the fa 
requirement. At low y, the hadronic system is close to the beam pipe and the fa' requirement 
affects the efficiency. The overall selection efficiency for CC events with Q 2 > 200 GeV 2 is 70%. 

The final sample consists of 1086 events. All events have been scanned visually, and no 
remaining cosmic or halo-muon background events have been found. The distribution of Q 2 
versus x for the accepted events is shown in Fig. 0. For x > 0.2, the sample is dominated by 
events for which the interaction vertex is determined by the calorimeter timing (full circles). 
Figure ^ also demonstrates that the acceptance is zero for events at very low x in the low Q 2 
region due to the fa cut. There is also zero acceptance for events at very low y (low Q 2 and high 
x), where a large part of the hadronic system escapes in the forward beam pipe direction. The 
MC is used to correct for the acceptance loss of such events in determining the cross sections. 

Figures H(a)-(d) show the distributions of the variables fa, 5, fa/Ex and 7^ in the final 
event sample, compared with the corresponding MC predictions, which include the contribu- 
tions from CC DIS and the small contributions from the background sources described below. 
The CC MC distribution is based on the CTEQ4D PDF set (see discussion in Sect. |9.1| ). 



The contamination from events not due to ep collisions, such as beam-gas interactions, is neg- 
ligible. In general, good agreement is observed, except for some excess at low 7^. Furthermore, 
the peak of the fa/Ex distribution is shifted to slightly lower values as compared to MC. For 
the events selected with a CTD vertex (Sect. |6^), shown by the open circles in Fig. ^ the 
distributions of Zvtx and the number of tracks assigned to the primary vertex are provided 
in Fig. |3](e) and (f). The good agreement between data and MC shows that as far as the 
acceptance calculation is concerned, the CC final state is well modeled by the MC. For the 
events selected with a timing vertex (Sect. |6~3|) , the distributions of Zyrx and the total energy 
in FCAL, Z?fcal, are shown in (g) and (h). About 18% of data events fall in this category, 
while MC predicts 16%. This small excess is directly related to the excess in the low j h (i.e. 
high x) region seen in (d). As seen in (h), all events in this sample have large FCAL energies 
to ensure a good resolution for Zvtx fr° m the timing. 

The relative resolution in Q 2 is approximately 20% over the entire range of Q 2 . The relative 
resolution in x improves from ~ 20% in the interval 0.01 < x < 0.0215 (see Sect. 0) to ~ 8% 
at high x, and that in y is approximately 8% over the entire range, except for y < 0.1 where 
it increases to ~ 11%. Here the resolutions are obtained by comparing the reconstructed 
quantities with the true values in MC, and the RMS value of the distribution is quoted. 

The fraction of background events in the final sample is typically below 1% at high Q 2 
and increases as Q 2 decreases, exceeding 10% in the interval 200 GeV 2 < Q 2 < 400 GeV 2 (see 
Sect. |7J), as estimated from MC. Photoproduction and Bethe-Heitler dilepton production in the 
dimuon channel are the dominant sources of background at low Q 2 , whereas the production 
and decay of time-like W bosons is the remaining background at high Q 2 . The contamination 
from NC events is negligible. 

7 Cross-section determination 

The single and double differential cross sections are determined using bin-by-bin unfolding. 
The measured cross section in a particular bin, cr mea s, is determined from 

-^obs ffibg (ir\\ 

C"meas J^C ' 



7 



where N ^ s is the number of observed events in the bin, iVb g is the estimated number of back- 
ground events, A is the acceptance and £ is the integrated luminosity. The acceptance, defined 
from the MC as the number of events reconstructed within the bin divided by the number of 
events generated in that bin, takes both the selection efficiency and the event migration due to 
resolution into account. 

The measured cross section includes the radiative effects discussed in Sect. |[ The correction 
factor to provide the Born-level cross section is defined as 



a 



SM 
Born 



Crad — SM ■ (11) 
a rad 



The numerator is obtained by numerically integrating ([|) over the bin. The value of af^, 
the cross section in the bin including radiation, is calculated using HERACLES 4.6.2 [3~I|. The 
measured Born-level cross section is then given by 

<^Born '^meas^-'rad • (1^) 

Finally, the quoted differential cross section, for example da/dQ 2 , is calculated as 



da 0"Born d(J SM 

dQ 2= ^L^Q 



Born 



(13) 



A similar procedure is used for da / dx , da / dy and d 2 a/dxdQ 2 . In this manner, the acceptance 
factor A is used to correct the effect of all the selection cuts (Sect. |) and the cross sections are 
extrapolated to the full kinematic range. In particular the MC is used to extrapolate beyond 
the ^/-region restricted by the yjB < 0.9 cut. The differential cross sections da/dx and da/dy 
are quoted in the region Q 2 > 200 GeV 2 . 

For the da/dQ 2 measurement, nine bins are used between Q 2 = 200 and 60000 GeV 2 . The 
bins have equal width in logioQ 2 between 400 and 22494 GeV 2 (four bins per decade), while 
the lowest and highest Q 2 bins have somewhat larger width. For da/dx, bins of equal width in 
logiox are used, three between x = 0.01 and x = 0.1, and four between x = 0.1 and x = 1.0. 
For the da/dy measurement, the y region is divided equally into two bins between y = 0.0 
and y = 0.2, and five bins between y = 0.2 and y = 0.9. All bins are defined such that their 
sizes significantly exceed the resolutions of the respective variables (see Sect. |T4|). The values 
at which da/dQ 2 and da/dx are quoted, Q 2 and x c , are chosen to be near the logarithmic 
center of each bin, except in the highest Q 2 and x bins, where they are chosen lower than the 
logarithmic center, reflecting the very steeply falling cross sections. The cross section da/dy is 
quoted at the center y c of each bin. 

The statistical errors are calculated using the square root of the number of measured events, 
N, for iV > 100 and otherwise from 68% Poisson confidence intervals around N. 

The values of Q 2 C1 x c and y c , the number of observed events, iV bs, the estimated number of 
background events, iVb g , the acceptance, A, and the radiative correction factor, C ra d, are given 
in Tables g|-@. 



8 Systematic uncertainties 

The major sources of systematic uncertainties for the quoted cross sections are described below. 
The individual uncertainties are added in quadrature separately for the positive and negative 
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deviations from the nominal cross section values to obtain the total systematic uncertainty. The 
uncertainty on the luminosity of the combined 1994-1997 sample is 1.6% and is not included 
in the total systematic uncertainty. 

• Uncertainty of the calorimeter energy scale 

The uncertainty in the hadronic energy scale of the calorimeter is determined by the 
methods described in to be 2% for FCAL and BCAL, and 3% for RCAL. Varying 
the energy scale of the calorimeter separately by this amount in the detector simulation 
induces shifts of the Jacquet-Blondel estimators for the kinematic variables. The resulting 
systematic errors in the measured cross sections are typically less than 10%, but increase 
to ~ 30% in the highest Q 2 bin and ~ 25% in the highest x bin. 

A 2% fraction of the accepted events have a measurable energy leakage from the CAL 
into the BAC. The average energy leakage for these events is 4 GeV. This effect on the 
cross-section measurement is negligible. 

• Variation of selection thresholds 

The threshold values of the selection cuts are varied independently in MC by typically 
10%. The largest effect is observed when varying the fa threshold; this changes the cross 
section by 8% in the lowest Q 2 bin and 4% in the highest y bin. Varying the fa' threshold 
produces a 3% change in the lowest y bin. 

• Uncertainty in the parton-shower scheme 

To test the sensitivity of the results to the details of the simulation of higher-order QCD 
effects in the hadronic final state, the lepto meps model is used instead of the ARIADNE 
model for calculating the acceptance. The largest effects are observed in the bins of lowest 
Q 2 (8%), highest Q 2 (6%), lowest x (5%) and highest y (5%). 

• Background subtraction 

The uncertainty of the photoproduction background is estimated by fitting a linear com- 
bination of the fa I Et distributions of the signal and the photoproduction MC samples to 
the corresponding distribution in the data, allowing the normalizations of the direct and 
resolved photoproduction components to vary. No cut on fa / E T is applied for this check. 
A 40% uncertainty in the photoproduction background is found, leading to a sizable sys- 
tematic error in the lowest Q 2 bins (4% in da/dQ 2 and maximum 7% in d 2 a/dxdQ 2 ). At 
Q 2 > 400 GeV 2 , the effect is less than 1%. 

• Trigger efficiency 

The simulation of the efficiency of the first-level trigger as a function of fa is examined 
by using data events triggered independently. This independent trigger is efficient for 
CC events and is based on the calorimeter energy sums. The difference between the 
efficiencies calculated from the data and from MC has a negligible effect on the measured 
cross section. 

• Choice of parton distribution functions 

The CC MC events are generated using the CTEQ4D PDFs. To examine the influence 
of variations of the PDFs on the cross-section measurement, the PDFs of MRS A |32|] and 
GRV94 |53 are also considered. In addition, a modification of the (i-quark to u-quark 



density ratio according to the prescription (d/u) + 5(d/u) = {d/u) + 0.1x(x + 1) (see |34[ 
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and discussion in Sect. 0) has been tested. This increases the predicted e + p CC cross 
section at high x. 

The MC events are re-weighted using these alternative PDFs and new acceptance correc- 
tion factors are computed. The change in the measured cross section is typically < 1% 
and at most 2.5% over the entire kinematic range of interest. 

The radiative correction factors also depend on the choice of the PDFs. Using the program 
hector f35fl , the difference in the radiative correction factors by changing the PDFs is 
found to be typically less than 0.1% and thus is neglected. 

The calculation of the differential cross section in a given bin according to (p~3|) uses the 
ratio of the differential to the integrated SM cross sections. The ratio is sensitive to the 
shape of the PDFs within this bin. The largest effects are observed in the highest Q 2 bin 
(3.5%) and highest x bin (3%). 

The effect of Fl 

The django program neglects the Fl contribution when generating CC events. The 
corresponding effect on the acceptance correction factors is evaluated by reweighting MC 
events with the ratio of the cross sections with and without Fl- The largest effect is 
observed in the highest y bin where it amounts to 1.5%. 

Uncertainty in the radiative correction 

The uncertainties in the radiative correction as determined with HERACLES are estimated 



to be smaller than 3% in the kinematic region considered f36| and are not included in the 
total systematic uncertainty. 



9 Cross-section results 

9.1 Single differential cross sections 

The differential cross sections da/dQ 2 , da/dx and da/dy are shown in Figs. f§(a), |5|(a) and ^|(a) 
and are compiledP] in Tables [3]-|3|. The uncertainty of the measured cross sections are dominated 
by statistical errors and are typically 7-20%. The cross section da/dQ 2 falls by a factor of 
about 50000 between Q 2 = 280 and 30000 GeV 2 . As a function of x, the cross section da/dx 
is largest at small x, showing a gradual decrease from x = 10 -2 to 2-10" 1 , followed by a sharp 
drop towards x = 0.6. The cross section da/dy decreases slowly as a function of y for y > 0.1. 

The Standard Model cross sections from (|2]) using CTEQ4D are also shown in Figs. f|-|], 
together with the ratios of the measured to the SM cross section. The Standard Model describes 
the data well, with the possible exception of da/dx, where the measurement at x > 0.3 is 
somewhat above the CTEQ4D prediction. 



Also shown are the cross-section predictions obtained from a NLO QCD fit [11], together 
with their uncertainties. The fit is made to the data from fixed-target experiments and NC 
DIS measurements at HERA in the region Q 2 < 5000 GeV 2 . Neither the CC data from HERA 
nor the recent high-Q 2 NC data from ZEUS || are included in the fit. The prediction from the 
fit also describes the data well, and in particular describes da/dx at higher x better than the 
CTEQ4D prediction^. The e + p CC DIS cross section is dominated by the d-quark contribution 

1 Tables || and ^ contain details of the systematic uncertainties that are correlated between cross-section bins. 
2 The CTEQ5 PDFs [|7|, which have been made available recently, predict higher da/dx than CTEQ4D in 
the region 0.1 < x < 0.6. 
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at high x, as can be seen from ([3]) and (|j) 



assumed has been of interest in recent years, for example see [34, 38 



The possibility of a larger d/u ratio than previously 

Modification || of PDFs 



with an additional term (S(d/u); see Sect. |8]) yields da/dx close to the NLO QCD fit as shown 
in Fig. |^(b). For comparison, the prediction of the MRST |3{J PDFs is also shown in Fig. [|(b). 



9.2 Double differential cross section 

The double differential cross section d 2 a/dxdQ 2 is compiled 1 in Table f|. The reduced double 
differential cross section, ex, is defined from the CC cross section as 



a 



G\ I M 2 W 



2ttx \ M 2 + Q 2 



-i 



' lWlL ' 1 [Y + F^(x,Q 2 )-Y^ c (x,Q 2 )-y 2 F c L c (x,Q 2 



dx dQ 2 2 

(14) 

The reduced cross sections as functions of x and Q 2 are displayed in Figs. |7| and [8|. The 
predictions of the CTEQ4D PDFs give a good description of the data, although at the highest 
x value the measured cross sections lie above the predictions. The predictions from the NLO 
QCD fit at high x are higher than those from CTEQ4D. 

In leading-order QCD, a depends on the quark momentum distributions as follows 

a = x [u + c+ (1 -y) 2 (d + s)] . (15) 

As a result, for fixed Q 2 , a at low x (i.e. high y) is mainly sensitive to the antiquark combination 
(u + c) while at high x (i.e. low y) it is dominated by the quark combination (d + s). These PDF 
combinations evaluated with the leading-order CTEQ4L parameterization are shown separately 
in Fig. [7|. Both the quark and the antiquark combinations are required in order to obtain a 
good description of the data. 



9.3 Comparison of NC and CC cross sections 

Fig. |9] compares the cross section da/dQ 2 for CC scattering with the ZEUS result for NC 
scattering ||. At low Q 2 , the CC cross section is much smaller than the NC cross section due 
to the relative strength of the weak force compared to the electromagnetic force. However, 
the CC cross section decreases with Q 2 less rapidly than that for NC scattering, reflecting the 
behavior of the W propagator as contrasted to the photon propagator which dominates NC 
scattering. At Q 2 ~ M^, M|, the CC and NC cross sections become comparable; at these 
large Q 2 , the weak and electromagnetic forces are of similar strengths. At yet higher Q 2 , the 
rapid fall of both CC and NC cross sections with Q 2 is due to the effects of the W and Z 
propagators, the decrease of the parton densities with increasing x, and, in particular for CC 
e + p scattering, the (1 — y) 2 term in the cross section. These observations were made in previous 
HERA measurements || [J ; here they are clearly demonstrated with a much improved precision. 



10 Electroweak analysis 

The absolute magnitude of the CC cross section, described by (fj), is determined by the Fermi 
constant Gf and the PDFs, while the Q 2 dependence of the CC cross section includes the 
propagator term [M^/(M^ + Q 2 )} 2 , which produces substantial damping of the cross section 
at high Q 2 . To compare the experimental results with the predictions of the Standard Model, 
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a x 2 fit to the measured differential cross section, da/dQ 2 , has been performed, treating Gp 
and M\y as free parameters purely for the purpose of this exercise. The result of this fit is 

G F = (l.l71±0.034(stat.)^32( s y st -)-aSi5(PDF)) x 10" 5 GeV" 2 , (16) 

and 

M w = 80.8ll| ( stat -)-43 (syst.)^;J (PDF) GeV. (17) 

The central values are obtained using the CTEQ4D PDFs. The major sources of systematic 
uncertainty in the determination of the cross section, namely the energy scale, the parton 
shower scheme and the luminosity, are taken into account in the systematic errors. The PDF 
errors quoted are obtained by re-evaluating the PDFs within the uncertainties given by the 
NLO QCD fit fTT |. The sensitivity of the result to variation of the value of Gp assumed in the 



extraction of the PDFs is negligible compared to the uncertainties quoted above. The point 
which gives minimum \ 2 (Xmm) * s displayed in Fig. [TO as the triangle together with the 70% 
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confidence level contour. This contour was determined using statistical errors only. The value 
of Gp obtained is in good agreement with the value Gp = (1.16639 ± 0.00001) x 10~ 5 GeV" 2 
obtained from muon decay ||10|| , implying the universality of the CC interaction over a wide 
range of Q 2 . The value of M w obtained agrees with the value of M w = 80.41 ± 0.10 GeV from 
the PDG fit [p7D| , using time-like production of W bosons at the Tevatron and at LEP. Since CC 



DIS represents space-like exchange, the result (|17]) is complementary to measurements of Mw 
from pp or e + e~ annihilation. This result constitutes an important experimental consistency 
check of the Standard Model. 

Two more fits are performed to determine Mw under more restrictive theoretical assump- 
tions. First, by evaluating the x 2 function along the line Gf = 1.16639 x 10" 5 GeV -2 , a mea- 
surement of the 'propagator mass' of the exchanged W boson can be made. The determination 
of the propagator mass is an important test of the SM description of the CC in the space-like 
regime. Second, a 'Standard Model fit' may be performed by evaluating the \ 2 function along 
the SM constraint 

M 2 1 
Vim-M^M^l-Ar' {Q) 

where Mz is the mass of the Z boson and a is the fine structure constant. The term Ar contains 
the radiative corrections to the lowest order expression for Gp and is a function of a and the 
masses of the fundamental bosons and fermions El . The constraint implied by flT8"D on Gp and 



Mw is also shown in Fig. [LC] as the heavy solid line. As can be seen from the plot, the value 
of Gf has a strong dependence on Mw- Therefore, within the context of the SM, the greatest 
sensitivity to Mw in this experiment may be obtained by combining the Mw dependence of 
the propagator term in the CC cross section with the M w dependence of Gp. 

The 'propagator-mass' fit to the measured differential cross section, da/dQ 2 , with Gp fixed 
to the value G F = 1.16639 x 10~ 5 GeV" 2 yields the result 

M w = 81.4l 2 ;^(stat.)±2.0(syst.)+3;o(PDF) GeV, (19) 



shown in Fig. [10] as the solid dot with horizontal error bars. The use of the Gf constraint has 
significantly reduced the uncertainty on the estimation of Mw compared to the value obtained 
in the unconstrained fit (|17]). The result is in agreement with the value of Mw obtained by 
direct measurement, M w = 80.41 ± 0.10 GeV fT~ 
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In order to use the SM constraint (|T8|), a, Mz, and all fermion masses, other than the mass 
of the top quark, M t , are set to the PDG values ||10|| . The central result of the fit was obtained 
with M t = 175 GeV and the mass of the Higgs boson Mh = 100 GeV. The x 2 function is 
evaluated along the line given by the SM constraint and the position of the minimum, shown 
in Fig. [10] as the large star, gives the 'Standard Model fit' result^ 



M w = 80.50lg:it(stat.)lS;S(syst.)±0.31(PDF)!JS(AM t , AM H , AM Z ) GeV. (20) 

The error labeled AM t , AMh, AMz is obtained by re-evaluating Mw with M t in the range 
170 < M t < 180 GeV, M H in the range 100 < M H < 220 GeV and M z in the range 91.180 < 
Mz < 91.194 GeV. The dependence of Mw on these changes is small, and the resulting error 
is negligible compared with the other errors quoted above. This result is in agreement with the 



value of Mw = 80.35 ± 0.21 GeV obtained in fixed-target neutrino-nucleon DIS ||41| . The good 
agreement with both direct and indirect determinations of Mw indicates that the SM gives a 
consistent description of a variety of phenomena over a wide range of energy scales. 



11 Summary 

Charged-current deep inelastic scattering, e + p —>■ v e X, has been measured with the ZEUS 
detector at HERA using 47.7 pb -1 of data collected during 1994-1997. Single differential 
cross sections da/dQ 2 , da/dx and da/dy have been presented with typical uncertainties of 
7-20%. The cross section da/dQ 2 falls by a factor of about 50000 as Q 2 increases from 280 
to 30000 GeV 2 . The double differential cross section d 2 a/dxdQ 2 has also been measured. 
A comparison between the data and Standard Model (SM) predictions shows clearly that 
contributions from antiquarks [u and c) and quarks (d and s) are both required by the data. 
The predictions of the SM, using recent parton distribution functions, give a good description 
of the full body of the data presented in this paper. 

The charged-current (CC) cross-section results for da/dQ 2 presented here have been com- 
pared with the recent ZEUS results for neutral-current (NC) scattering. The CC cross section 
is found to fall with Q 2 less rapidly and to approach the NC cross section at Q 2 > M w , M|, in 
agreement with previous observations at HERA. This shows that the weak and electromagnetic 
forces reach similar strengths for Q 2 above M Wl M\. 

A fit to the data for da/dQ 2 with the Fermi constant Gf and Mw as independent parameters 
yields 

G F = (1.171±0.034 (stat.)+J;^ (syst.)+^ (PDF)) x 10" 5 GeV' 2 , 
M w = 80.8+^5 (stat.)5° (syst.)^ (PDF) GeV. 

A propagator fit, with Gp fixed to the PDG value yields 



M w = 81.4+ 2 ;£(stat.)±2.0(syst.)+3;o(PDF) GeV. 

These results show that the SM gives a consistent description of charged-current induced time- 
like and space-like processes over a wide range of virtualities. 



3 It should be clearly understood that ( |20| ) represents a constrained fit using assumptions on the validity of 
the Standard Model which implies that the result cannot be used in a global average of experimental values of 
M w . 
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Table 1: The differential cross section da/dQ 2 for the reaction e + p — > v e X. The following 
quantities are given for each bin: the Q 2 range; the value at which the cross section is quoted, 
Q 2 ; the number of selected events, iV bs; the number of expected background events, iVb g ; the 
acceptance, A; the radiative correction factor, C rad (see Sect. ^); the measured Born-level cross 
section, da/dQ 2 ; and the Born-level cross section predicted by the Standard Model, using 
CTEQ4D PDFs. The first error of each measured cross section value gives the statistical error, 
the second the systematic uncertainty. 



Q 2 range 


Ql 


N obs 


N bg 


A 


Crad 


da/dQ 2 (pb/GeV 2 ) 


(GeV 2 ) 


(GeV 2 ) 


measured 


SM 


200 - 400 


280 


141 


16.5 


0.46 


1.01 


2.94 ±0.28 ^ -10~ 2 


2.80 


io- 


-2 


400 - 711 


530 


173 


5.3 


0.64 


1.01 


1.82 ±0.14 ±0.08 -10- 2 


1.87 


io- 


-2 


711 - 1265 


950 


248 


2.3 


0.74 


1.01 


1.29 ±0.08 ±0.03 -10- 2 


1.15 


io- 


-2 


1265 - 2249 


1700 


205 


2.3 


0.79 


1.03 


5.62 ±0.40 ±0.08 -10^ 3 


6.07 


io- 


-3 


2249 - 4000 


3000 


169 


2.4 


0.81 


1.04 


2.62 ±0.20 1^9 T0~ 3 


2.61 


io- 


-3 


4000 - 7113 


5300 


91 


0.6 


0.82 


1.05 


7 Q1 +0.93 +0.38 1fl -4 
'• Ji -0.83 -0.31 ' iu 


8.29 


io- 


-4 


7113 - 12649 


9500 


45 


0.7 


0.86 


1.07 


2.00 ±0.17 -10" 4 


1.65 


io- 


-4 


12649 - 22494 
22494 - 60000 


17000 
30000 


13 
1 


0.1 
0.01 


0.96 
1.48 


1.07 
1.09 


o ci +0^95 +0.45 i n -5 
z - ul -0.72 -0.38 ' ±u 
cr q +14. +1.8 in-7 
°- J -4.9 -1.5 ±U 


1.71 
6.24 


io- 

10" 


-5 
-7 
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Table 2: The differential cross section da/dx for the reaction e + p — > z/ e X for Q 2 > 200 GeV 2 . 
The following quantities are given for each bin: the x range; the value at which the cross 
section is quoted, x c ; the number of selected events, N ohs ; the number of expected background 
events, N hg ; the acceptance, A; the radiative correction factor, C rad (see Sect. 0); the measured 
Born-level cross section, da/dx; and the Born-level cross section predicted by the Standard 
Model, using CTEQ4D PDFs. The first error of each measured cross section value gives the 
statistical error, the second the systematic uncertainty. 



x range 





N obs 


Nb 9 


A 


C-rad 


0.015 


136 


4.0 


0.56 


1.04 


0.032 


246 


7.4 


0.79 


1.02 


0.068 


306 


10.0 


0.85 


1.02 


0.130 


200 


2.9 


0.82 


1.01 


0.240 


124 


0.6 


0.73 


1.01 


0.420 


46 


0.06 


0.56 


1.00 


0.650 


3 


0.00 


0.36 


1.02 



da/dx(ph) 



measured 



±0^0™ 
±0.17±0.07 
±0.08 ±0.03 
±0.49i^ 2 7 
±0.23 US 
tli ±0.6 

+7.8 +2.1 
-4.4 -1.6 



SM 



0.0100 - 0.0215 
0.0215 - 0.0464 
0.0464 - 0.1000 
0.1000 - 0.1780 
0.1780 - 0.3160 
0.3160 - 0.5620 
0.5620 - 1.0000 



4.50 

2.64 

1.44 

6.88 

2.57 

6.8 

8.1 



10 2 


3.97 


10 2 


10 2 


2.76 


10 2 


10 2 


1.50 


10 2 


10 1 


7.04 


10 1 


10 1 


2.39 


10 1 


10° 


4.56 


10° 


io- 1 


3.55 


io- 1 
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Table 3: The differential cross section da/dy for the reaction e + p — ► v e X for Q 2 > 200 GeV 2 . 
The following quantities are given for each bin: the y range; the value at which the cross section 
is quoted, y c \ the number of selected events, N ohs ; the number of expected background events, 
N hg ; the acceptance, A; the radiative correction factor, C rad (see Sect. the measured Born- 
level cross section, da/dy; and the Born-level cross section predicted by the Standard Model, 
using CTEQ4D PDFs. The first error of each measured cross section value gives the statistical 
error, the second the systematic uncertainty. 



y range 


Vc 


N obs 


Nb 9 


A 


Crad 


0.00 - 0.10 


0.05 


192 


7.3 


0.64 


0.98 


0.10 - 0.20 


0.15 


249 


5.2 


0.85 


0.99 


0.20 - 0.34 


0.27 


240 


2.6 


0.83 


1.01 


0.34 - 0.48 


0.41 


185 


5.0 


0.78 


1.03 


0.48 - 0.62 


0.55 


117 


3.8 


0.70 


1.05 


0.62 - 0.76 


0.69 


61 


2.1 


0.60 


1.07 


0.76 - 0.90 


0.83 


42 


1.1 


0.44 


1.10 



da/dy(pb) 



measured 



SM 



6.95 ±0.52 ±0.28 
5.92 ±0.38 ±0.14 
4.27 ±0.28 toil 
3.52±0.27±0.08 
2.46 ±0.24ig-{jjj 



i cp;+0.23 
I. OO _ a20 



1.49 



+0.27 
-0.24 



+0.07 
-0.10 

±0.13 



•10 1 
•10 1 
•10 1 
•10 1 
•10 1 
•10 1 
•10 1 



6.92 
5.79 
4.30 
3.11 
2.34 
1.84 
1.53 



10 1 
10 1 
10 1 
10 1 
10 1 
10 1 
10 1 
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Table 4: The double differential cross section da/dxdQ 2 for the reaction e + p — > v e X. The following quantities are given for each bin: 
the x and Q 2 range; the values at which the cross section is quoted, x c and Q 2 C ; the number of selected events, N ohs ; the number of 
expected background events, N hg ; the acceptance, A; the radiative correction factor, C rad (see Sect. 7); the measured Born-level cross 
section, da/dxdQ 2 ; and the Born-level cross section predicted by the Standard Model, using CTEQ4D PDFs. The first error of each 
measured cross section value gives the statistical error, the second the systematic uncertainty. 



Q 2 


range 




x range 


n 2 


x c 


N obs 


N bg 


A 


C 




d 2 a/dxdQ 2 (pb/GeV 2 ) 


(GeV 2 ) 








l bev ) 


rad 


measured 


SM 


onn 
zUU 


— 4UU 


n 


ni nn 
U1UU 


n noi k 
— U.UZlO 


98n 

ZoU 


U.UIO 


ou 


9 n 
z.u 


n rzl 

U.O'l 


i 
i 


nn 

UU 


7 n 
I .U 


-1.0 ±u - 7 


i n-i 
•1U 


K 07 
O.Z ( • 


1 n- 
1U 


-1 






n 
u 


noi k 


n n/ifi/i 
— U.U404 




n n^9 
u.uoz 


49 


A 7 


n sn 

U.oU 


n 
u 


QQ 

yy 


1 nn 
l.yy 


+0.40 +0.19 
-0.34 -0.18 


1 n-i 
•11) 


Z.U i ■ 


i n- 


-1 






n 
u 


U404 


n 1 nnn 
— U.1UUU 




n nfis 

U.UOo 


98 
zo 


Q 

o.y 


n 7i 

U. 1 1 


n 
u 


Q7 

y i 


fi o 

o.z 


+1.8 i-i n 
-1.5 ±Lu 


i n — 2 
•1U 


7 a n 
/.4y • 


1 n- 
1U 


-2 


4DD 


- 71 1 

1 J. J. 


n 
u 


ni nn 
UlUU 


n noi ^ 
— U.UZlO 


c^n 


n ni ^ 

U.UIO 


AO 


1 

1.0 


n p.a 


i 
i 


n9 
uz 


3 86 


+0.72 +0.33 
-0.62 -0.39 


•10- 1 

1U 


4 1 n • 


i n~ 


-1 






n 
u 




- D n464 




n n^9 
u.uoz 


^9 
oz 


U.o 


n 87 

U.o ( 


n 
u 


QQ 

yy 


1.62 


±0.09 


•io- 1 


1.76 • 


10" 


-1 






o 


0464 


- 0.1000 




0.068 


44 


1.0 


0.86 





98 


6.4 


+1.1 +0.3 
-1.0 -0.2 


•io- 2 


6.60- 


10" 


-2 









1000 


- 0.1780 




0.130 


21 


0.5 


0.77 





98 


2.41 


+0.67 +0.13 
-0.53 -0.07 


•10~ 2 


2.44- 


10" 


-2 


711 


- 1265 





0100 


- 0.0215 


950 


0.015 


37 


0.3 


0.50 


1 


06 


2.97 


+0.58 +0.30 
-0.49 -0.25 


■io- 1 


2.72 • 


10" 


-1 









0215 


- 0.0464 




0.032 


85 


0.4 


0.88 


1 


01 


1.53 


+0.19 +0.04 
-0.17 -0.07 


•io- 1 


1.32 • 


10" 


-1 









0464 


- 0.1000 




0.068 


69 


0.5 


0.89 





99 


5.59 


+0.76 +0.12 
-0.68 -0.11 


•IO" 2 


5.31 • 


io- 


-2 









1000 


- 0.1780 




0.130 


43 


0.2 


0.84 





98 


2.61 


+0.46 , q in 
-0.40 ^U-J-U 


•IO" 2 


2.02- 


10" 


-2 









1780 


- 0.3160 




0.240 


14 


0.0 


0.66 





99 


5.7 


Hi ±0.3 


•IO" 3 


5.90- 


10" 


-3 


1265 


- 2249 





0215 


- 0.0464 


1700 


0.032 


58 


0.9 


0.77 


1 


05 


6.9 


+1.1 +0.3 
-0.9 -0.2 


•IO" 2 


7.90- 


10" 


-2 









0464 


- 0.1000 




0.068 


66 


1.1 


0.91 


1 


02 


3.03 


+0.43 +0.05 
-0.38 -0.07 


•IO" 2 


3.63- 


10" 


-2 









1000 


- 0.1780 




0.130 


40 


0.3 


0.89 


1 


00 


1.31 


+0.24 +0.01 
-0.21 -0.03 


•io- 2 


1.47- 


10" 


-2 









1780 


- 0.3160 




0.240 


28 


0.3 


0.83 


1 


00 


5.1 


+L2 ±0 1 
-1.0 ^u- 1 


•io~ 3 


4.42- 


10" 


-3 



Table 



4: continued. 



Q 2 range 


x range 


Ql 




N obs 




A 


Crad 


d 2 a/dxdQ 2 (pb/GeV 2 ) 


(GeV 2 ) 






(GeV 2 ) 


measured 


SM 


2249 - 4000 


0.0464 - 


0.1000 


3000 


0.068 


76 


1.1 


0.88 


1.06 


2 


i 9 +0.28 +0.02 
—0.25 —0.11 


•10- 2 


1.92 


10" 


-2 




0.1000 - 


0.1780 




0.130 


42 


0.2 


0.89 


1.03 


8 


q +1.4 +0.2 
— 1.2 —0.1 


•10~ 3 


8.86 


10" 


-3 




0.1780 - 


0.3160 




0.240 


27 


0.1 


0.86 


1.01 


2 


7q +0.65 +0.08 
' —0.54 —0.09 


•10- 3 


2.84 


10" 


-3 




0.3160 - 


0.5620 




0.420 


15 


0.0 


0.75 


1.00 


9 


r +3.2 +0.4 
° —2.4 —0.6 


•10~ 4 


5.03 


10" 


-4 


4000 - 7113 


0.0464 - 


0.1000 


5300 


0.068 


22 


0.4 


0.69 


1.07 


6 


3 t\i ±0.4 


•10- 3 


7.65 


10" 


-3 




0.1000 - 


0.1780 




0.130 


35 


0.2 


0.92 


1.04 


3 


ye +0:75 +0.23 
' -0.63 -0.19 


•10- 3 


3.86 


10" 


-3 




0.1780 - 


0.3160 




0.240 


26 


0.1 


0.89 


1.03 


1 


r-i +0.36 +0.07 
01 -0.30 -0.09 


•10- 3 


1.43 


10" 


-3 




0.3160 - 


0.5620 




0.420 


8 


0.0 


0.84 


1.04 


2 


7 +1.3 +0.3 
' -1.0 -0.2 


•10- 4 


2.73 


10" 


-4 


7113 - 12649 


0.1000 - 


0.1780 


9500 


0.130 


15 


0.6 


0.84 


1.08 


1 


06 t° f 8 ±0.13 


•10- 3 


1.05 


10" 


-3 




0.1780 - 


0.3160 




0.240 


18 


0.1 


0.94 


1.05 


5 


6 til ±0- 4 


•10- 4 


4.74 


10" 


-4 




0.3160 - 


0.5620 




0.420 


10 


0.0 


0.90 


1.05 


1 


82 ±g"P ±0.19 


•10- 4 


1.09 


10" 


-4 


12649 - 22494 


0.1780 - 


0.3160 


17000 


0.240 


5 


0.1 


1.04 


1.07 


8 


n +5:5 +1.7 
u -3.5 -1.5 


•10" 5 


7.97 


10" 


-5 




0.3160 - 


0.5620 




0.420 


7 


0.0 


1.03 


1.06 


6 


+3.4 +1.0 
-2.3 -0.9 


•10~ 5 


2.72 


10" 


-5 



Table 5: The differential cross sections da/dQ 2 , da/dx and da/dy for the reaction e + p — > v e X. 
The following quantities are given for each bin: the value at which the cross section is quoted; 
the measured Born-level cross section; the statistical uncertainty; the total systematic uncer- 
tainty; the uncorrelated systematic uncertainty and those systematic uncertainties with sig- 
nificant (assumed 100%) correlations between cross-seciton bins. The systematic uncertainties 
considered to be correlated were: the FCAL energy scale (Si); the BCAL energy scale (S2) and 
the uncertainty in the parton-shower scheme (63). 



da/dQ' 2 


Ql 


da/dQ 2 (pb/GeV 2 ) 


<W (%) 


5 sys (%) 


Knc (%) 


Si 


(%) 


52 


(%) 


53 


(%) 


280 
530 
950 
1700 
3000 
5300 
9500 
17000 
30000 


2.94 • 10- 2 
1.82 • 10- 2 
1.29- 10- 2 
5.62 • 10~ 3 
2.62 • 10~ 3 
7.91 • 10- 4 
2.00 • 10~ 4 
2.61 • 10~ 5 
5.9 • 10- 7 


±9.5 

±7.8 

±6.4 

±7.1 

±7.8 
+12. 
-11. 

+18. 

-15. 

+37. 

-28. 
+233. 
-84. 


±12. 

+4.2 
-4.5 
+2.3 
-2.2 
+1.5 
-1.4 
+1.3 
-3.3 
+4.8 
-3.9 

±8.6 

+17. 
-15. 
+30. 
-26. 


+8.4 
-8.3 
+3.6 
-3.7 
+2.0 
-1.8 
+1.0 
-1.1 
+0.8 
-1.2 
+0.8 
-0.6 

±0.7 

+1.3 
-1.0 
+3.2 
-3.3 


-1.0 
+2.4 
-1.4 
+1.0 
+0.3 
+0.5 
+0.4 
-0.2 
+0.6 
-2.1 
+3.0 
-1.7 
+4.4 
-4.5 
+7.1 
-5.6 
+9.0 
-9.4 


-1.0 
+1.9 
-1.2 
+0.7 
-1.0 
+1.0 
+0.5 
+0.0 
+0.4 
-2.2 
+3.6 
-3.4 
+7.1 
-7.0 
+15. 
-13. 
+28. 
-23. 


-7.9 
+7.9 
-1.9 
+1.9 
-0.5 
+0.5 
+0.9 
-0.9 
-0.8 
+0.8 
+0.4 
-0.4 
+2.2 
-2.2 
+3.3 
-3.3 
+7.6 
-7.6 


da/dx 


x c 


da/dx (pb) 


tfstat (%) 


Ssys (%) 


Sunc (%) 


Si 


(%) 


S2 


(%) 


S3 


(%) 


0.015 
0.032 
0.068 
0.130 
0.240 
0.420 
0.650 


4.50 • 10 2 
2.64- 10 2 
1.44- 10 2 
6.88- 10 1 
2.57- 10 1 
6.8-10° 
8.1 • 10- 1 


±8.8 
±6.6 
±5.9 
±7.2 
±9.0 

+17. 
-15. 
+97. 
-54. 


+7.6 
-7.9 
+2.7 
-2.6 
+1.9 
-2.2 
+2.4 
-1.8 
+3.5 
-3.7 
+8.3 
-8.4 
+26. 
-20. 


±6.1 

+2.4 
-2.5 

±1.5 
±1.2 
±1.0 

±2.4 

+5.2 
-5.1 


-1.0 
+1.5 
-0.2 
+0.2 
-0.4 
+0.5 
+0.4 
-0.3 
+2.8 
-3.2 
+7.8 
-7.9 

+25. 
-19. 


-2.3 
+0.6 
-0.3 
+0.7 
+0.4 
-1.2 
+2.0 
-1.3 
+1.8 
-1.6 
+0.9 
-1.1 
+2.1 
-0.9 


-4.3 
+4.3 
-0.7 
+0.7 
-1.0 
+1.0 
-0.1 
+0.1 
+0.2 
-0.2 
+1.5 
-1.5 
+1.3 
-1.3 


da/dy 


Vc 


da/dy (pb) 


<Wt (%) 


S sys (%) 


<5unc (%) 


Si 


(%) 


s 2 


(%) 


S3 


(%) 


0.05 


6.95 • 10 1 


±7.5 


±4.1 




±3.7 




-0.8 
+0.9 




-0.5 
+0.1 




-1.3 
+1.3 


0.15 
0.27 
0.41 
0.55 
0.69 
0.83 


5.92- 10 1 
4.27- 10 1 
3.52 • 10 1 
2.46 • 10 1 
1.55 • 10 1 
1.49 • 10 1 


±6.5 
±6.5 
±7.6 
±9.6 

+15. 
-13. 
+18. 
-16. 


±2.4 

+2.5 
-2.2 
+2.2 
-2.4 
+3.6 
-3.2 
+4.6 
-6.7 
+8.9 

-8.8 




±1.2 

+1.9 
-1.8 
+2.1 
-2.2 
+2.8 
-2.6 
+3.2 
-3.4 
+5.6 
-5.8 




+0.3 
-0.5 
+1.2 
-0.8 
+0.5 
-0.2 
+0.4 
+0.2 
-0.9 
-1.0 
+1.7 
+0.5 




-0.8 
+0.8 
-0.5 
+0.9 
+0.3 
-0.6 
+2.0 
-1.4 
+3.1 
-5.6 
+5.8 
-5.7 




-1.8 
+1.8 
-0.9 
+0.9 
-0.2 
+0.2 
-1.1 
+1.1 
-1.0 
+1.0 
-3.4 
+3.4 
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Table 6: The double differential cross section da/dxdQ 2 for the reaction e + p — > v e X. The following quantities are given for each bin: 
the Q 2 and x values at which the cross section is quoted, Q 2 C and x c ; the measured Born-level cross section; the statistical uncertainty; 
the total systematic uncertainty; the uncorrelated systematic uncertainty and those systematic uncertainties with significant (assumed 
100%) correlations between cross-seciton bins. The systematic uncertainties considered to be correlated were: the FCAL energy scale 
(5i); the BCAL energy scale (5 2 ) and the uncertainty in the parton-shower scheme (5 3 ). 



Vcl LreV ) 


x c 


(1 O 1 CIXLU^ ^po/vorev ) 


fstat I/O ) 


S ys I/O ) 


<Junc I/O J 


0\ I/O J 


2 l/oj 




ZoU 


n m k 

U.U10 


( .u ■ 1U 


+17. 

-15. 


+9.2 
-9.9 


ztO.O 


-1.8 
+1.3 


-3.6 
+1.1 


-7.2 
+7.2 


ZoU 


U.UoZ 


1 qq 1 n — 1 
i.yy ■ ru 


+20. 
-17. 


+9.6 
-9.1 


+3.5 
-3.9 


-1.1 
+2.3 


+0.6 
+3.0 


-8.1 
+8.1 


ZoU 


U.UOo 


c o i n— 2 
O./ • 1U 


+29. 
-24. 


4-1 7 
±11. 


±10. 


-1.4 
+3.1 


+0.3 
-0.1 


-5.0 
+5.0 


oou 


n m ^ 

U.U10 


O.oO • 1U 


+19. 

-16. 


+8.6 
-10. 


+7.9 
-8.4 


-2.3 
+0.8 


-3.2 
-1.0 


-3.4 
+3.4 


CQfl 

OoU 


U.Uoz 


l.uz • ru 


+16. 
-14. 


1 r A 
±0.4 


±D.U 


-0.4 
-1.0 


-1.6 
+1.9 


-0.6 
+0.6 


OOKJ 




6 4-1 fr 2 

U . L ± 1U 


+18. 
-15. 


+4.5 
-3.7 


+1.9 
-2.2 


-1.9 
+3.2 


+0.5 
+0.8 


-2.3 
+2.3 


530 


0.130 


2.41 ■ 10~ 2 


+28. 
-22. 


+5.5 
-3.1 


+3.8 
-2.3 


-0.4 
+2.9 


+1.5 
+0.5 


+1.9 
-1.9 


950 


0.015 


2.97- lO" 1 


+19. 
-17. 


+10. 
-8.5 


+8.8 
-8.2 


+1.7 
+3.5 


-1.6 
+2.8 


-1.4 
+1.4 


950 


0.032 


1.53 • 10 1 


+12. 
-11. 


+2.8 
-4.4 


+1.9 
-2.0 


-0.4 
+0.0 


-3.3 
+0.4 


+2.0 
-2.0 


950 


0.068 


5.59- 10~ 2 


+14. 
-12. 


+2.2 
-2.0 


+1.6 
-1.0 


-0.1 
-0.1 


-1.0 
+0.6 


-1.4 
+1.4 


950 


0.130 


2.61 • 10~ 2 


+18. 

-15. 


+3.8 
-3.6 


+1.4 
-1.3 


+0.6 
+0.2 


-0.8 
+1.2 


-3.3 
+3.3 


950 


0.240 


5.7- 10~ 3 


+35. 
-27. 


+5.4 
-4.9 


+3.2 
-3.1 


+1.3 
+0.0 


+1.5 
+0.4 


-3.9 
+3.9 


1700 


0.032 


6.9- 10~ 2 


+15. 
-13. 


+4.3 
-3.1 


+3.2 
-2.6 


-0.1 
+0.8 


+2.1 
+0.9 


+1.8 
-1.8 


1700 


0.068 


3.03- 10~ 2 


+14. 
-13. 


+1.6 
-2.2 


+0.5 
-0.8 


+0.6 
+0.2 


-1.6 
+0.0 


+1.4 
-1.4 


1700 


0.130 


1.31 • 10~ 2 


+19. 

-16. 


+1.1 
-1.9 


+0.8 
-1.5 


-0.2 
-0.7 


+0.8 
-1.0 


+0.2 
-0.2 


1700 


0.240 


5.1 • 10" 3 


+23. 
-19. 


+2.6 
-1.9 


+1.1 
-0.7 


+2.1 
-1.7 


+1.1 
-0.1 


+0.5 
-0.5 



Table 6: continued. 



Q 2 (GeV 2 ) 


x c 


d' 2 a/dxdQ 2 (pb/GeV 2 ) 


<Wt (%) 


S sys (%) 


Sunc (%) 


Si (%) 


s 2 (%) 


h (%) 


3000 


0.068 


2.12 


• 10~ 2 


+13. 

1 


+1.1 

— O.l 


+0.6 
— i.y 


-1.3 
— i.y 


+0.5 

A 1 
— 4. 1 


-0.8 

-hU.o 


3000 


0.130 


8.0 • 


10~ 3 


+18. 
— io. 


+1.9 

— 1.0 


+0.6 
— U.O 


+1.1 

— L.O 


+0.7 

+ 1.U 


-0.8 

-hU.o 


3000 


0.240 


2.79 


• 10~ 3 


+23. 
— 19. 


+2.8 
—3.2 


+0.9 
—0.8 


+2.5 
—2.9 


-0.8 
+0.6 


+0.3 
—0.3 


3000 


0.420 


9.5 • 


10 -4 


+33. 

— ZD. 


+4.5 

— O.O 


+0.9 

1 A 


+3.9 

— O. i 


+0.4 
— u.y 


-2.1 

1 O 1 

+Z.1 


5300 


0.068 


6.3 • 


10~ 3 


+27. 
—22. 


+5.5 
—5.8 


+1.6 
— 1.3 


+1.5 
+2.2 


+4.3 
—5.4 


-1.5 
+1.5 


5300 


0.130 


3.75 


■ 10~ 3 


+20. 

1 1 
— 1 1 . 


+6.1 

— 0.1 


+1.3 

n a 
— U.l 


+1.3 
— U.o 


+5.2 

A A 


+2.5 
— z.o 


5300 


0.240 


1.51 


• 10~ 3 


+24. 
—20 


+4.9 
—5 8 


+0.6 

—0 7 


+4.8 
—5 7 


+0.8 
— 1 i 


-0.3 
+0 3 


5300 


0.420 


2.7- 


10" 4 


+49. 
-35. 


+9.8 
-7.6 


+1.3 
-0.5 


+9.6 
-7.6 


+0.6 
+0.9 


+0.0 
+0.0 


9500 


0.130 


1.06 


• io- 3 


+35. 
-27. 


±12. 


+1.4 
-1.5 


+0.7 
-1.3 


+12. 
-11. 


+3.8 
-3.8 


9500 


0.240 


5.6- 


io~ 4 


+30. 
-23. 


+6.4 
-7.7 


±0.7 


+4.9 
-5.4 


+3.7 
-5.2 


+1.6 
-1.6 


9500 


0.420 


1.82 


• io~ 4 


+43. 
-31. 


±11. 


+0.5 
-0.6 


±11. 


-0.3 
-0.2 


-0.5 
+0.5 


17000 


0.240 


8.0 • 


IO" 5 


+69. 
-44. 


+21. 
-19. 


+1.8 
-1.7 


+1.0 
-1.2 


+20. 
-19. 


+3.5 
-3.5 


17000 


0.420 


6.3 • 


IO" 5 


+54. 
-37. 


+16. 
-14. 


±0.9 


+14. 
-13. 


+4.6 

-3.3 


+3.5 
-3.5 




Figure 1: (a) A schematic diagram of charged-current positron-proton scattering, (b) A view of 
a charged-current candidate event in the ZEUS detector, projected in the plane parallel to the 
beam. The filled boxes indicate energy deposits in the calorimeter. The transverse momentum 
imbalance can be clearly seen in the calorimeter and also from the tracks of charged particles 
measured in the central tracking detector. 
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ZEUS CC 1994-97 




Figure 2: Distribution of the selected CC candidates in the (x, Q 2 ) plane. Open (full) circles 
show the events selected with (without) tracking vertex. The curve shows the rfr cut of 12 GeV. 
The bin boundaries are shown by the dotted lines, delimited by the diagonal dotted line of the 
kinematic limit, y — 1. The bins used in the double differential cross section measurement are 
marked with solid lines. 
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ZEUS CC 1994-97 




Figure 3: Comparison of the final CC data sample (solid points) with the expectations of 
the MC (histograms), normalized to the luminosity of the data. The distributions of (a) the 
missing transverse momentum, Ifr, (b) the variable 5, (c) Ifx / Et, the ratio of missing transverse 
momentum to total transverse energy and (d) the variable 7^, are shown. In (e) and (f), the 
distributions of the Z position of the event vertex and the number of tracks assigned to the 
primary vertex, Nyr, are shown for selected events with CTD vertex (see Sect. |6.2| ). In (g) and 
(h), the distribution of Zvtx and -Efcal are shown for events passing the selection with the 
timing vertex (see Sect. pOl ). 
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ZEUS CC 1994-97 
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Figure 4: (a) The e + p CC DIS Born cross section da/dQ 2 for data (solid points) and the 
Standard Model (SM) expectation evaluated using the CTEQ4D PDFs. (b) The ratio of the 
measured cross section da/dQ 2 to the SM expectation evaluated using the CTEQ4D PDFs. 
The statistical errors are indicated by the inner error bars (delimited by horizontal lines), 
while the full error bars show the total error obtained by adding the statistical and systematic 
contributions in quadrature. Also shown by a dot-dashed line is the result of the NLO QCD 
fit together with the associated PDF uncertainties (shaded band). 
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ZEUS CC 1994-97 




Figure 5: (a) The e + p CC DIS Born cross section da/dx for data (solid points) and the 
Standard Model (SM) expectation evaluated using the CTEQ4D PDFs. (b) The ratio of the 
measured cross section da/dx to the SM expectation evaluated using the CTEQ4D PDFs. 
The statistical errors are indicated by the inner error bars (delimited by horizontal lines), 
while the full error bars show the total error obtained by adding the statistical and systematic 
contributions in quadrature. Also shown by a dot-dashed line is the result of the NLO QCD 
fit together with the associated PDF uncertainties (shaded band). The dashed line represents 
the result of modifying the d/u ratio with S(d/u) = 0.1x(x + 1). The dotted line shows the 
prediction from MRST PDFs. 
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ZEUS CC 1994-97 




Figure 6: (a) The e + p CC DIS Born cross section da/dy for data (solid points) and the 
Standard Model (SM) expectation evaluated using the CTEQ4D PDFs. (b) The ratio of the 
measured cross section da/dy to the SM expectation evaluated using the CTEQ4D PDFs. 
The statistical errors are indicated by the inner error bars (delimited by horizontal lines), 
while the full error bars show the total error obtained by adding the statistical and systematic 
contributions in quadrature. Also shown by a dot-dashed line is the result of the NLO QCD 
fit together with the associated PDF uncertainties (shaded band). 
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ZEUS CC 1994-97 




Figure 7: The reduced cross section b as a function of x, for fixed values of Q 2 . The dots 
represent the data, while the expectations of the Standard Model evaluated using the CTEQ4D 
PDFs are shown as the solid lines. For illustration, the leading-order PDF combinations x(u+c) 
and (l—y) 2 x(d+s), taken from the CTEQ4L parameterization, are plotted as dotted and dashed 
lines, respectively. Also shown is the result of the NLO QCD fit (dash-dotted line). 
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ZEUS CC 1994-97 
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Figure 8: The reduced cross section a as a function of Q 2 , for fixed values of x. The dots 
represent the data, while the expectations of the Standard Model evaluated using the CTEQ4D 
PDFs are shown as the solid lines. Also shown is the result of the NLO QCD fit (dash-dotted 
line) . 
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ZEUS e + p DIS 1994 - 97 




CT (GeV*) 



Figure 9: Comparison of the differential cross sections da/dQ 2 for NC (solid dots) and CC(open 
dots) deep inelastic e + p scattering from the ZEUS 1994-97 analysis. The lines represent the SM 
predictions evaluated using CTEQ4D PDFs. The second highest-Q 2 bin in the NC measurement 
contains no event: thus an upper limit of the cross section at 95% confidence level is quoted 
and indicated by the arrow head. This cross section is quoted at Q 2 = 30400 GeV 2 (very close 
to the CC point at Q 2 = 30000 GeV 2 ), as indicated by the lower end of the bent line. 
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Figure 10: The result of a fit of the CC DIS cross section to determine Gf and My/. The 
triangle indicates the result of the fit, i.e. the position of Xmm> the minimum value of \ 2 ■ The 
70% confidence level contour is shown as the ellipse. The dot with error bar shows the result 
of the 'propagator-mass' fit, described in the text, in which the \ 2 function is evaluated along 
the horizontal line Gf = 1.16639 x 10" 5 GeV" 2 . The SM constraint implied by fll8|) is shown 
as the heavy solid line. The large star shows the position of Xmm(SM), the minimum of the x 2 
function evaluated along the SM constraint line. The solid bars crossing the SM constraint line 
show where x 2 (SM) = ^(SM) + 1. The small star at the bottom with the error bar shows 
the value of My/ obtained in the 'Standard Model fit'. Note that all errors and the confidence 
level contour correspond to statistical errors only. Also shown (dotted) are lines of constant 
G F (1.16639 x 10" 5 GeV" 2 ) and constant M w (80.41 GeV). 
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